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THE ENERGETIC ENVIRONMENT AND THE DENSE INTERSTELLAR 
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RESUMEN 

Presentamos los ajustes a la distribution espectral de energfa, desde el cercano-infrarrojo hasta radio, de 30 
galaxias luminosas y ultraluminosas en el infrarrojo con modelos que incluyen formacion estelar (SB) y actividad 
nuclear (AGN). La motivation de este trabajo fue determinar importantes parametros fisicos para esta clase de 
objetos, como son: la profundidad optica hacia la fuentc infrarroja, el ritmo y la eficiencia de la formacion estelar 
y la contribution del AGN. Encontramos que, aunque la mitad de la muestra necesita modelos combinados 
(i.e. AGN+SB), solo el 30% necesito un AGN contribuyendo mas del 10% a la luminosidad infrarroja. En 
todos los casos la contribution del SB cs dominantc. De nucstros modelos tambicn pudimos obtcner la masa de 
gas molecular denso. Si asumimos que esta masa es la misma que la que cs trazada por la molccula dc HCN, 
podemos rcproducir la relation lineal cntre la luminosidad dc HCN y la masa dc gas denso observada por Gao 
y Solomon (2004). Sin embargo, el factor de conversion entrc la luminosidad dc HCN y la masa de gas denso 
que derivamos es un factor 2 menor que el asumido por estos autores. Por ultimo, tambien encontramos que 
la eficiencia de la formacion estelar disminuyc a mcdida que el brote de formacion estelar envejece. 

ABSTRACT 

We fit the near-infrared to radio spectral energy distributions of a sample of 30 luminous and ultra-luminous 
infrared galaxies with models that include both starburst and AGN components. The aim of the work was 
to determine important physical parameters for this kind of objects such as the optical depth towards the 
luminosity source, the star formation rate, the star formation efficiency and the AGN fraction. We found that 
although about half of our sample have best-fit models that include an AGN component, only 30% have an 
AGN which accounts for more than 10% of the infrared luminosity whereas all have an energetically dominant 
starburst. Our models also determine the mass of dense molecular gas. Assuming that this mass is that 
traced by the HCN molecule, we reproduce the observed linear relation between HCN luminosity and infrared 
luminosity found by Gao & Solomon (2004). However, our derived conversion factor between HCN luminosity 
and the mass of dense molecular gas is a factor of 2 smaller than that assumed by these authors. Finally, we 
find that the star formation efficiency falls as the starburst ages. 

Key Words: Ulirgs — ISM: dust extinction — Galaxies: active — Infrared: galaxies 



1. INTRODUCTION 

With total infrared luminosities between 10 11 — 
10 12 L/0 and > 10 12 L , respectively, Luminous and 
Ultraluminous Infrared galaxies, (U)LIRGs, are the 
most luminous objects in the local universe. Al- 
though scarce at low redshift, they may account for 
the bulk of all star formation activity at z > 2 — 3 and 
dominate the far-infrared background (e.g. Blain 
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et al. 2002). Many of them are found in merg- 
ing systems (e.g. Sanders ct al. 1988), suggesting 
that dynamical interaction has driven gas towards 
the nucleus, fueling a massive starburst (SB) and/or 
the central massive black hole (Mihos & Hernquist 
1996). Despite extensive investigation over the last 
decades, there is still considerable uncertainty as to 
the nature of their power source. In order to analyze 
the power mechanism in these sources, tracers that 
do not suffer large extinction have been preferred 
recently, such as hard X-rays, and MIR to radio di- 
agnostics. The weakness of the hard X-ray luminos- 
ity seen in most ULIRGs have been interpreted by 
Risaliti et al. (2006) as a real lack of AGN activ- 
ity. MIR diagnostics such as the strengths of PAH 
emissions, the 9.7/zm absorption feature, the high 
ionization lines, and the MIR continuum slope indi- 
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cate that the starburst is dominant in about 80% of 
ULIRGs. However, there are strong discrepancies in 
the results obtained with the different methods (e.g. 
Armus et al. 2007). This is in line with the recent 
finding that the strength of PAHs and the shape of 
the mid-infrared continuum cannot be safely used to 
disentangle AGN and SB contributions (Vega et al. 
2005, Weedman et al. 2005). In the radio domain, 
the most direct way of distinguishing between AGN 
and starburst power sources in ULIRGs is to search 
for very compact radio continuum emission towards 
the nuclei. Nagar et al. (2003), by using 15 GHz 
radio continuum data with a resolution of 150 mas, 
concluded that most of the 83 ULIRGs of their sam- 
ple are AGN powered. However, Smith et al. (1998) 
detected radio-SNe in the NW nucleus of Arp 220 in 
a region of "2 x "A and concluded that no AGN is 
necessary to explain the IR luminosity in this source. 

We followed a different approach to examine the 
energetic environment in (U)LIRGs. Instead of con- 
sidering a single spectral region, we base our study 
on the thorough analysis of the panchromatic spec- 
tral energy distribution (SED). Thus, we first deter- 
mine the observed NIR - radio SED of a sample of 30 
local (U)LIRGs and then, these data are compared 
with suitable models with SB and AGN components. 

A Similar approach has been adopted by Farrah 
et al. (2003). They used SB and AGN templates 
to model the SED from the optical-submm range. 
We prefer to use models built on well calibrated star 
forming complexes instead of templates re-scaled to 
different luminosities because it avoids the dubious 
process of "re-scaling" a "template" to another lu- 
minosity, and because the use of models gives access 
to a more sound physical picture of the environment 
within which such objects evolve, that can be used 
for the study of even more extreme phases such as 
those likely found at high redshift (e.g. Mortier et 
al. 2005), which is one of the ultimate goals of the 
present investigation. We excluded from the analysis 
the optical data because, due to the large extinction, 
they are likely dominated by the underlying old stel- 
lar population and not by the starburst (e.g. Farrah 
et al. 2001). In contrast, the inclusion of radio data 
in the analysis gives us an independent constraint on 
the SB strength and age (Bressan, Silva & Granato 
2002). 

2. SAMPLE 

The selection criteria for our sample of 30 
(U)LIRGs are summarized as follows: (a) The galax- 
ies are compact (U)LIRGs from Condon et al. (1991) 
(i.e. brighter than 5.25 Jy at 60 [m\) with radio data 
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at 3 or more frequencies, (b) We also require the 
galaxies to have a well-sampled infrared SED, with 
data in the NIR, MIR and FIR spectral ranges. For 
all the sources of our sample we collected NIR data 
from the 2MASS All-Sky Extended Source Catalog 
and IRAS data from IRS A and NED, respectively 
Radio data were taken from Clemens et al. (2008). 
Additional IR and sub-mm broad-band data were 
taken from the literature. We also found Spitzer 
IRS spectra in the SSC archive for eleven objects 
of our sample. Our objects are compact and isolated 
so that we are confident that our fluxes sample the 
entire starburst in all cases. 

3. MODELS 

In order to find the best fit model we adopt the 
following approach: we first test the null-hypothesis 
that our sample sources are pure SB. Then, if the ob- 
served SED shows an excess emission over the best 
fit in those wavelength regions that are more affected 
by the presence of an embedded AGN, typically in 
the NIR-FIR range, we reject the null-hypothesis 
and add an AGN component to the model. The 
AGN model is selected from a library generated with 
the radiative transfer code developed by Granato & 
Danese (1994) and added to the SB model in such a 
way that the total luminosity is 

L{ ot = (l-/)xLf +/xLt GN ,with0 < / < 1. (1) 

Notice that in this work our aim is just to assess the 
possible contribution of the warm AGN component 
to the total SED of ULIRGs, without any pretense 
of discriminating the precise AGN properties. The 
best fit is obtained by minimizing the merit function, 
X, calculated as 

1 ( F mod{i) ~ F obs (i) \ 2 

x = n2- ( { EMi) ) ' (2) 

where F mo d(i), F b s (i), and Err(i) are the model 
flux values, the observed fluxes and observational er- 
rors respectively. N is the number of passbands used 
for the fit. 
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Fig. 1. Examples of the comparison between the broad band SED (open circles) and the best fit model (thick solid line) 
for six galaxies of our sample. The additional emission from the AGN, when needed, is indicated by a dot-dashed line. 
When available, the IRS low resolution spectrum (solid thick line) and the corresponding SED model (solid thin line) 
are shown upward displaced by one dex. The residuals between models and data are shown in the bottom panels as 

(./"data /"model)/ /data- 



The SB models are selected from a library gener- 
ated with GRASHl] (Silva et al. 1998) and involve 

6 GRASIL can be run via the WEB interface 
web.pd.astro.it/galsynth/ which is maintained at OAPD, 
INAF by L. Paoletti, A. Petrella & D. Selvestrel, 
and can also be downloaded from the webpage 
http: / /adlibitum. oats. inaf.it/silva/default. html 



seven main free parameters: the e-folding time and 
age of the burst, t\, and age^, the sub- mm dust spec- 
tral index, /?, the escape time, t esc , the molecular 
cloud (MC) optical depth at 1 fim, n, the fraction 
of gas in MCs, f mc , and the core radius of the King 
profile for the dust-star distributions, r. The range 



4 



VEGA, CLEMENS, BRESSAN, ET AL. 



of values of these free parameters are listed in Table 
[1] For a complete description of models we refer the 
reader to Vega et al. (2008) and references therein. 

4. RESULTS AND CONLUSIONS 

In Fig. [T]wc show the SED fits for six galaxies of 
our sample. The corresponding best-fit parameters 
of the SB components for all the sample are shown in 
Table [2] The last four columns of the table display 
the value of the evolutionary phase of the SB, Esb, 
the current star forming rate, SFR, the mass of the 
dense component, M^cm and the fractional contri- 
bution of the AGN to the 8-1000 /jm flux indicated 
as AGN/Tot. The SED fits for all the sample can be 
seen in Vega ct al. (2008)0. In Fig[2]wc plot all the 
SED best-fit models of our sample of (U)LIRGs nor- 
malized to the 60 /jm flux density. The plot shows 
the large heterogeneity of (U)LIRGs SEDs. In par- 
ticular, Arp 220, which is usually used as prototype 
of ULIRGs, has one of the most extreme SED. 
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the PAH equivalent widths and mid-infrared spec- 
tral slope overestimate (often by large factors) the 
AGN fraction. Therefore, our work supports the 
notion that neither the mid-infrared slope nor the 
PAH equivalent widths provide a good estimate of 
the AGN contribution (see conclusions in Vega et al. 
2005). For two galaxies, our findings are markedly 
different to those of Armus et al., namely IR 14348- 
1447 and IR 15250+3609. For the former, we find 
an AGN contribution of ~ 10%, while they find no 
evidence of an AGN from the MIR analysis. In- 
deed, this AGN contribution is required by its red 
NIR colour, J-K = 1.71 which cannot be explained 
by stellar populations alone. The power source of 
IRAS 15250+3609 could not be determined by Ar- 
mus ct al., while our analysis indicates that it is 
a mixed object with a very obscured AGN, that 
contributes ~ 20% of the IR luminosity. This is 
where the combined analysis of the IRS spectra and 
panchromatic SED is most powerful. 
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Fig. 2. Comparison of SED models for our sample of 
(U)LIRGs normalized to the 60 fmi flux density. 

4.1. Star formation versus AGN 

The 53% of the sources show evidence for the 
presence of an AGN. However, only in 9/30 sources 
the contribution of the AGN to the total IR luminos- 
ity exceed 10%. Only one galaxy, IR 08572+3915, 
seems to be completely powered by an AGN. We 
therefore find that the fraction of SB dominated ob- 
jects in our sample (i.e. pure SB galaxies and mixed 
galaxies with AGN contributions < 10%) is 70%. 
Our measured AGN fractions are consistent with 
those reported by Armus ct al. (2007) based on the 
MIR line ratios Ne[V]/Nc[II] and 0[IV]/Ne[II]. But, 

7 SED models and data are available in electronic form at 
http:/ /www. aanda.org and they can also be downloaded from 
http:// www. adlibitum. oats . inaf . it/ silva/def ault .html 



Fig. 3. Correlation between the dense gas mass, derived 
from our models, and the AGN corrected IR luminosity 
for the 30 (U)LIRGs. The thick solid line is the linear 
fit. The dotted line corresponds to the fit between the 
IR luminosity and the dense molecular mass as traced by 
HCN emission obtained by GS04. The cross in the upper 
left denotes the typical errors of the IR luminosities and 
masses derived with our models 

4.2. Molecular masses 

The MIR-FIR spectral region in our (U)LIRGs is 
dominated by the emission from the molecular clouds 
where the star formation is taking place. Therefore 
we can determine quite accurately the mass in dust 
associated with the star forming regions via the ex- 
tinction required to reproduce the infrared part of 
the SED. The dust mass can be converted into molec- 
ular gas mass by using a typical dust/gas ratio (i.e. 
G/D =100). We find a relation between the infrared 
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Galaxies with an * have Spitzer spectra (projects IDs 105, and 14, P.I.: J. R. Houck, and ID 30323, P.I.: L. Armus). 

4.3. The star formation efficiency 



luminosity and molecular gas mass with the same 
slope as that derived by Gao & Solomon (2004, here- 
after GS04) from observations of HCN emission, sec 
Fig El However, we derived a constant of conversion 
between the HCN luminosity and the mass of dense 
molecular gas a factor of 2 smaller than that assumed 
by these authors for their global sample of galaxies, 
but in agreement with the conversion factor found 
by Gracia-Carpio et al. (2006) for their (U)LIRG 
population. Actually, an accurate determination of 
the conversion factor between HCN and the mass of 
the very dense molecular gas needs more extensive 
studies, including the combined analysis of the SED 
and the high excitation HCN transitions (Vega et al. 
2009 in prep.). 



With the ratio between the SFR and the Mden, 
we obtain a quantity that measures the reciprocal 
of the SF time scale and that is usually called star 
formation efficiency (SFE, GS04). We do not find 
clear correlations between SFE and IR luminosity of 
the SB, or the mass of molecular gas. However we 
do find an evident anti-correlation between SFE and 
the SB phase (^sb = 2 f s S see Vega et al. 2005, 
for a thorough discussion on the SB phases). The 
SFE decreases by about one order of magnitude as 
the SB evolves from its early phase to the more ad- 
vanced phases (upper panel of Fig. 0]). Such a rapid 
decrease of the SFE as the SB evolves suggests that 
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Fig. 4. Upper panel: SFE as derived from our models ver- 
sus the evolutive phase of the starburst. Galaxies in com- 
mon with the GS04 sample are marked as stars. Filled 
symbols correspond to galaxies fitted with pure SB mod- 
els while open symbols correspond to those fitted with 
mixed models. Lower panel: Correlation between the 
far- infrared color (/60//100) and the derived SFE for our 
sample of (U)LIRGs. Crosses denote the pure SB galax- 
ies, while open circles denote the galaxies with AGN. The 
filled circles denote the mixed galaxies but with the IR 
colours corrected for the AGN contributions. 

feedback plays a major role in the evolution of these 
objects. £"sb is a quantity derived from models and 
it would be more desirable to put the above finding 
on observational grounds. Vega et al. (2005) found 
a tight correlation between E§b and the FIR colour 
/60//100 i n the sense that as the SB evolves, i.e. E$b 
increases, the /60//100 decreases. Thus, by using the 
FIR colour instead of E$b , we obtain the correlation 
shown in the lower panel of Fig. [4J 

\og(SFE) = (1.74 ± 0.2) log(-p-) + (1.55 ± 0.03) 

(3) 

GS04 claim that the SFE is constant for all galaxies 
and that the star formation law in terms of dense 



molecular content is well represented by a power law 
with exponent ~1.0. Instead of SFE we could have 
used the ratio Lir/Lhcn and obtained a correlation 
analogous to that shown in Figure 9a of GS04. Con- 
sidering only ULIRGs in that figure, a definite trend 
of decreasing efficiency at cooler /60//100 colour, 
with slope ~ 1.72, is seen. Therefore, this corre- 
lation between the observed colour and the observed 
SFE indicates that the variation of the SFE is real. 
The additional information on the evolutive phase 
of the SB and on the AGN contamination, allowed 
by a realistic modelling procedure, eliminates a large 
fraction of the dispersion in the SFE, found by GS04. 
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